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CREATING  EQUATION  OF  STATE  OF 
CONDENSED  SYSTEMS  FROM  DYNAMIC 
EXPERIMENTS 

V.  Ye.  Fortov 
Moscow 

Successive  theoretical  calculation  of  the  thermodynamic 
properties  of  & substance  in  the  condensed  state  involves  the 
quantum  mechanics  problem  of  many  bodies  when  there  exists  a 
strong  interaction  among  them.  At  the  present  time  this  cannot 
be  achieved  in  a general  form.  In  this  situation  experimental, 
primarily  dynamic,  research  methods  based  on  obtaining  the 
necessary  state  by  compression  and  irreversible  heating  of  the 
medium  at  the  front  of  a strong  shook  wave  [lj  acquire  particular 
significance . 

The  characteristic  feature  of  dynamic  experiments  is  direct 
registration  of  only  the  mechanical  parameters  of  shock  compres- 
sion - pressure  (p),  specific  volume  (V),  and  internal  energy  (E). 
To  construct  the  complete  equation  of  state  from  these  data  we 
must  introduce  model  representations  of  the  properties  of  the 
studied  substance.  Here  the  data  of  the  experiment  were  used 
to  determine  the  constants  and  functions  contained  in  the  equation 
of  state;  for  certain  terms  purely  theoretical  expressions  are 
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often  used.  The  criterion  of  correctness  for  equations  of 
state  thus  obtained  is  agreement  between  calculated  and  measured 
states  of  shock  compression.  Here  a situation  may  arise  in  which 
various  equations  of  state  result  in  different  values  of  the 
Gruneisen  coefficient  and  temperature  for  the  shock-compressed 
substance,  even  though  the  original  experimental  material  was 
well  described  [2]. 

The  thermodynamically  complete  equation  of  state,  however, 
may  be  plotted  from  dynamic  experiments  without  introducing  limit- 
ing ideas  concerning  the  properties  and  the  nature  of  the  studied 
medium  [3,  4]. 

In  dynamic  experiments  registration  of  the  kinematic  param- 
eters of  propagation  of  a shock  discontinuity  makes  it  possible, 
by  using  the  general  laws  of  conservation  on  the  shock  wave 
front,  to  determine  the  equation  of  state  in  a caloric  form 
I>E(p,  V)  for  the  entire  region  of  the  pV  plane  accessible  in 
the  experiment  [1].  The  internal  energy,  however,  does  not 
represent  thermodynamic  potential  with  respect  to  the  pV  variable, 
and  in  order  to  plot  the  complete  thermodynamics  of  the  studied 
system  we  need  the  additional  dependence  T*T(p,  V). 

iT’-'Co  ?u  i nr  from  the  first  principle  of  IhcrmodytMnil  »■;;  and 
assuming  that  from  the  experiment  the  dependence  i>E(p,  V)  Is 
known,  we  car;  easily  obtain  [3] 

(i> 

The  solution  to  this  linear  heterogeneous  differential 
e . nation  with  variable  coefficients  is  the  dependence  T»T(p,  V) 
which  interests  us.  Equation  (1)  has  a corresponding  system  of 
characteristic  differential  equations: 
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the  first  of  which  determines  pressure,  the  second  - temperature 
on  the  isentrope.  Equations  (1)  or  (2)  and  (?)  are  completed  by 
boundary  conditions:  temperature  must  be  given  in  the  region 

where  it  can  bo  experimentally  determined.  In  the  present  cal- 
culations the  temperature  was  assigned  as  a function  of  specific 
volume  on  the  isobar  of  atmospheric  pressure;  the  experimental 
data  uf  v«V(T)  with  respect  to  the  thermal  expansion  of  the 
substance  were  used. 

The  caloric  equation  of  state  E«E(p,  V)  was  constructed  from 

N 

experimental  points  , pit  chaiotically  arranged  in  the 

pV  plane  in  the  form  of 


*<**)•  %£'***+  (U) 

In  (^)  the  coefficients  are  found  by  the  method  of  least 
squares  [4].  In  determining  coefficients  e^  the  weight  inversely 
proportional  to  registration  error  of  the  shock  compression 
states  was  assumed  by  each  experimental  point. 

Th"  degree  of  the  polynomial  q in  (4)  was  selected  by 
analysing  the  residual  sum  of  squares  when  the  degree  Increased 
q-*q+l.  Here  the  statistical  Fischer  significance  criterion  [5] 
was  used. 

Thv  equation  of  state  in  caloric  form  (^)  can  be  used 
directly  for  hydrodynamic  calculations  of  adiabatic  flows. 
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isentropes  are  determined  by  Integrating  differential  equations 
-)  and  (3^;  The  shock  adiabats  are  found  by  numerical  solution 

* too  Rar.kin-Kugono  1 1 relationship: 

y(p+/v)(iv  \)-n*(n.  V)  F.c(p*.V,).  (*>) 

• V)  is  taken  in  the  form  of  (4). 

The  dependence  tetween  the  accuracy  of  the  obtained  solution 
a -rro;  in  th.  original  data  was  determined  by  the  Monte  Carlo 
hod  by  simulating  the  probability  structure  of  the  measuring 
..cess  cn  the  electronic  computer  [4],  The  computer  reassembled 
ssifcle  combinations  of  random  factors  resulting  in  experimental 
nor  and  determined  the  effect  of  these  factors  on  the  solution 
. quat ions  (2  ■ and  (3).  The  proposed  method  enabled  a compre- 
oslve  estimate  for  the  accuracy  of  the  solution  obtained  for 
riven  configuration  of  experimental  points  in  the  pV  plane 
u specific  level  of  experimental  errors. 


•a;  an  cxi;;.in.-  we  Introduce  the  results  for  the  equation  of 
ate  for  nir/.cl.  For  this  we  used  experimental  data  for  dynamic 
r:  cession  of  .’.olid  [6,  ?,  8]  anu  porous  specimens  of  a given 
.ement . In  odnitica  the  internal  energy  (according  to  [10])  was 
•:s : gnel  at  tv.  a points  on  the  atmospheric  pressure  isobar.  As 
original  uata  for  system  (2),  (3)  we  use  the  results  of 
• siring  ‘.he  thermal  expansion  of  nickel  [11]  at  p»l  bar 


.Sortie  a:r.  lysis  [^]  showed  th.*1:  to  describe  the  existing 
• or:  inental  data  -..a  must  select  the  degree  of  the  polynomial 
. ; c=3.  Figure  1 shows  the  Hugonoit  adiabats  (Ni)  plotted 
.■■or.;  the  obtained  equation  of  state.  Tne  difference  between 
to  cl:  compression  states  calculated  according  to  (U)  and  ih) 
vd  those  which  were  measured  lies  within  the  limits  of  experimental 
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Figure  1.  HugonoiL  ad.' abate  for 
nickel  of  different  initial  poros- 

PH 

ttiea  m» — (p  - density  of  crystal, 
p0 

p - density  of  original  specimen); 

p - zero  isentrope  for  equations 

(2)  and  (*0;  dashed  line  - zero 
isentrope  for  [2];  0— [ 7 J ; V—  C 8 3 ; 

A~[6J;  §-p»l  bar). 

An  important  characteristic  of  the  substance  used  in  all 
model  theories  for  the  equation  of  state  is  the  "cold"  compression 
curve  PX(V)  - the  dependence  of  pressure  on  specific  volume  at 
T»0°K.  Existing  methods  for  obtaining  this  curve  from  the  results 
of  the  experiment  involve  introducing  a large  number  of  theoretical 
ideas  about  the  properties  of  the  studied  substance.  Assuming 
that  the  isotherm  T*0°K  coincides  with  the  isentrope,  the  cold 
compression  curve  in  this  approach  can  be  obtained  by  Integrating 
equation  (2)  for  TQ*0°K.  The  elastic  pressure  curve  thus  plotted 
is  shown  in  Pig.  1,  where  it  is  compared  with  px(V)  from  [2], 

It  should  bo  mentioned  that,  in  contrast  to  existing  theories, 
there  was  no  preliminary  division  of  the  thermodynamic  factors 
into  "thermal"  and  "cold"  components  in  the  method  of  constructing 
an  equation  of  state  which  is  discussed  here. 

Temperature  was  calculated  by  combined  integration  of 
equations  (2)  and  (3)*  For  the  initial  value  of  the  specific 
volume  temperature  was  determined  by  means  of  the  thermal 
expansion  curve  introduced  into  the  computer  (p»l  bar).  With 
these  initial  conditions  system  (2),  (3)  was  integrated  up  to 
the  Intersection  of  the  isentrope  with  the  shock  adlabat  in  the 
pV  plane  of  the  solid  substance.  Figure  2 shows  the  results 
of  determining  temperature;  here  we  also  see  for  comparison  the 
results  of  calculating  temperatures  on  shock  adiabats  for  solid 
nickel  according  to  the  Mie-Ortineisen  theory  [6]  and  in  the 


approximation  of  the  theory  of  free  volume  [ 7 U - 

The  characteristic  temper- 
ature "step"  on  the  shock 
adiabat  is  the  result  of  the 
melting  of  the  substance  at 
the  front  of  the  shock  wave. 

The  absence  of  detailed  ex- 
perimental data  in  the  region 
of  thte  melting  curve  resulted 
in  a situation  in  which  the 
melting  curve  coincided  with 
the  isentrope  in  this  rase; 
the  real  temperature  "stop" 
(according  to  [12])  is  located 
somewnere  above  that  which  is  shown  in  Fig.  2. 

The  statistical  analysis  made  by  Monte  Carlo  method  showed 
that  if  accuracy  in  the  nature  of  the  arrangement  of  the  experi- 
mental points  in  the  pV  plane  are  considered  the  error  in  deter- 
mining temperature  on  the  Hugonoit  adiabat  for  solid  nickel  is 
.).n  the  order  of  'j-8% , which  is  independently  confirmed  by 
comparison  to  temperature  calculations  for  model  theories  of  the 
b q ■ a i* . o n s t a t e 
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Figure  2.  Calculation  for 
temperature  of  nickel:  1 - from 

equations  .(2),  (3),  C21);  2 - 
study  [6  j ; 3 - study  [7];  s - 
1 sent ropes . 

Key:  (a)  Liquid;  (b)  Liquid* sol- 

id; (c)  Solid. 
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